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Q^', Based on a sample of approximately 500,000 hadronic Z° decays accumulated 

between 1993 and 1998, the SLD experiment has set limits on 24 fully charged 
two-body and quasi two-body exclusive charmless hadronic decays of B^, B^, 
and B^ mesons. The precise tracking capabilities of the SLD detector provided 
for the efficient reduction of combinatoric backgrounds, yielding the most precise 



^ ' available limits for ten of these modes. 
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The search for exclusive charmless decays of B mesons is motivated by the CKM- 
suppression of the VT-boson mediated h ^ u transition, which suppresses the leading order 
weak decay to charmless final states by a factor of |KtbP/|KfeP — 10~^ relative to that of 
charmed final states. Thus, observation of exclusive charmless modes with even modest 
branching fractions can indicate the participation of heretofore unobserved physical pro- 
cesses. 

Recently, several results have increased the interest in exclusive charmless B meson 
decays. The CLEO collaboration |jl| has improved their measurement of the decay B — ^ 
TjK*, with the measured branching fractions Br{B^ -^ rjK*^) = (2.73:^o|2 ± 0.50) x 10~^ 
and 5r(5° -^ r]K*^) = (1.38l|]i^ ± 0.17) x 10~^ somewhat above the expected range 
of (0.02 - 0.82) X 10"^ and (0.01 - 0.89) x 10-^ respectively. In addition, the DELPHI 
collaboration has reported a measurement [Q] of the combined mode Br{B^ —* p^ix^ + 
K^'^-K^) = (ITl^g^ ± 2) X 10~^, again somewhat higher than both the expected range of 
(0.4 — 2.0) X 10^^ and the corresponding CLEO measurements |Q of Br{B^ -^ p^tt"*") = 
{l.5toi ± 0.4) X 10-5 and Br{B+ -^ K*%+) < 2.7 x 10"^ 

In this Letter, we present limits from the SLD detector on several two, three, and four- 
prong fully charged two-body and quasi-two-body final states. Although the B meson sample 
available to the SLD detector is fairly limited in comparison to those produced at LEP, 
CESR, and the TEVATRON, the excellent tracking and a priori knowledge of the B meson 
production point admit limits competitive with those produced elsewhere. Most limits 
presented here on four-prong final states, for which combinatoric backgrounds are worst, are 
the first available. In addition, the cms energy available to experiments running at the Z° 
pole allows the study of Bs decays, which are inaccessible to experiments running at the 
T(4^). 

Search modes reported in this Letter include B^, 5° -^ P^P^ (two-prong), B^ -^ P^V'^ 
(three-prong), and B^^B^^ -^ V^V^ (four-prong), and their charge-conjugates, where P = 
n,K is a. stable pseudoscalar meson and V = p°, K*^, is a vector meson resonance with a 
sizeable branching fraction into two charged pseudoscalar mesons (100%, 66.7%, and (49.1 
± 0.8)% p, respectively). The ability to fully reconstruct the decaying B meson state, with 
precise momentum and vertex information for each of the charged daughter tracks, provides 
an essential constraint in the analysis; no attempt was made to search for modes with one 
or more long-lived final state neutral particles. 

The SLD detector instruments the sole interaction region of the SLAC Linear Collider 
(SLC). The luminous region of the SLC is an ellipsoid of dimensions approximately 2 and 
0.8 fim in the horizontal (x) and vertical (y) directions perpendicular to the beam axis, and 
700 nm along the beam axis. Due to motion of the collision point, however, the location of 
the luminous region is known to only ~ 7 fim in x and y. 

Charged particle tracks are reconstructed in the central drift chamber (CDC) and the 
CCD-based pixel vertex detector (VXD) in a uniform axial magnetic field of 0.6T. Including 
the uncertainty in the location of the luminous region (IP) of the SLC, the VXD2 vertex 
detector, in place through 1995, exhibited an r — (r — z) impact parameter resolution of 
11/xm (38/im) at high momentum, and 71pm (80pm) at p^V^nO = 1.0 GeV/c 0. The 
corresponding resolution for the VXDS vertex detector J3, in place since 1996, is 14pm 
(26.5pm) at high momentum, and 33pm (33pm) at p±y/sm9 = 1.0 GeV/c. The combined 
CDC/ VXD momentum resolution in the plane perpendicular to the beam axis is Sp±/p± = 



.01)2 _l_ (^.0026p_L/GeV/c)2. High momentum charged tracks are reconstructed in the range 
|cos6'| < 0.85, with an efficiency of 96% for |cos6'| < 0.65. A segmented Si-W forward 
calorimeter, with polar angle acceptance between 23 mr and 68 mr, is used to monitor the 
SLC luminosity via t-channel Bhabha scattering. 

The SLD accumulated an integrated luminosity of 19.1 pb^^ of e~^e~ annihilation data at 
the Z° pole between 1993 and 1998. Of this, 14.0 pb^^ was taken with the upgraded VXD3 
vertex detector in place. 

The complete reconstruction of the fully-charged final state provides a number of con- 
straints which can be used to discriminate between signal and potential background sources. 
Candidate track combinations must be consistent with having arisen from a single vertex. 
This vertex is displaced from the collision point by an average of ~3mm, which is measured 
with an average uncertainty of 60/im for the search modes. The point of closest approach 
of the extrapolated vertex momentum resultant to the SLD IP ('vertex impact parameter') 
must be consistent with zero. The invariant mass of the tracks forming the vertex must 
be consistent with that of the B meson and have a total momentum consistent with known 
fragmentation properties. Tracks emerging from the 5-meson decay vertex should have a 
relatively small opening angle, a large momentum, and a relatively large impact parameter 
with respect to the SLD IP. For quasi two-body modes involving vector meson resonances 
{B —>■ PV and B —>■ VV), two of the charged tracks must have an invariant mass consistent 
with that of each resonance. In addition, for B -^ PV modes, the decay angle 6h of the 

V state with respect to its flight direction ('hehcity angle') must be consistent with the 
distribution dT/dQh oc cos^ 6h dictated by angular momentum conservation. 

Candidate decays were reconstructed by considering all combinations of two, three, and 
four tracks which pass track quality cuts and with total charge for PP and VV candi- 
dates, and ±1 for PV candidates. The invariant mass of the candidate decay was required 
to be above 5.05 GeV/c^ (5.15 GeV/c^) for B+ and 5° (5°) modes. The probability of the 
vertex fit to the candidate tracks was required to be greater than 1.0% (0.5% for B — > VV 
modes and B -^ PV with V = (p), with a significance (separation from the SLD IP divided 
by the associated error) of greater than 1.0 (0.6 for B — * VV modes and B —>■ PV with 

V = (p). The smallest impact parameter D, normalized to its corresponding uncertainty, of 
any track in the candidate vertex was required to be greater than 1.1 (0.6 for B -^ VV modes 
and B -^ PV modes with V = (p). The change in the vertex invariant mass between the 
assumption of the nominal (kaon) mass and pion mass for all relevant tracks in the candidate 
vertex was required to be less than 0.3 GeV/c^ (1.2 GeV/c^) for the B -^ PK^ {B -^ Kp°) 
modes; this cut suppresses background vertices which get an artificially large mass due to 
a mistaken mass hypothesis for one or more tracks. A second mass reconstruction quantity 
J^vv: defined to be the sum of the absolute values of the differences between the recon- 
structed and nominal masses of the vector meson and B meson candidates, was required to 
be less than_0.6 GeV/c^, 0.4 GeV/c^, and 0.4 GeV/c^ for the B -^ p°p°, B -^ p^V {V ^ pO), 
and B^ — > K*^K*^modes, respectively. For relevant modes, the reconstructed vector meson 
masses were required to be in the ranges [0.2 — 1.1], [0.7 — 1.0], and [1.000 — 1.035] GeV/c^ 
for V = p°, K*^, and 0. Finally, | cos6h\ was required to be greater than 0.3 for B -^ PV 
modes. 

To further suppress background, an ad-hoc discriminator function was devised, and tuned 
to a sample of Monte Carlo (MC) Z^ -^ bb events approximately ten times that of data, and 



a sample of light quark (udsc) events approximately four times that of data. For B -^ PP 
modes, this function took the form 
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with m the invariant mass of the candidate vertex, S the vertex significance, A the largest 
angle between any tracks belonging to the vertex (A = 1/A — 0.9), P the vertex fit probability, 
D the minimum normalized impact parameter, / the vertex impact parameter (/ x 1000 + 
24 = 1/(J + 0.001)), and X = E^ert/Ebeam the scaled vertex energy. AM is the difference 
in the vertex mass between the pion and kaon hypotheses, and exploits the propensity for 
all tracks deriving from decays of the various search modes to be at high momentum. The 
parameters Cj > 0, mg, /q, and Iq were tuned separately for the individual search modes, 
while Mb was set to 5.28 GeV/c^ for S° or B+ and 5.37 GeV/c^ for B^. 
For B -^ PV modes, the discriminator function took the form 



TpV = J^Vp\a,a-n + ^^^ ^^'"'"^' + (1 " COs(/l7r)), 



with m„ the invariant mass of the vector meson candidate, and h = cos Oh- The vector meson 
masses My were set to 0.77, 0.89 and 1.02 GeV/c^ for p, K*^ and candidates, respectively, 
with corresponding widths 6m^ of 0.1, 0.08 and 0.006 GeV/c^, respectively. 
For B -^ VV modes, the discriminator function took the form 
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with vector meson candidates selected according to the track partition yielding vector meson 
masses closest to those of the search mode. 

The discriminator functions were tuned for the various search modes by maximizing the 
separation between the SLD MC sample (which contains no charmless hadronic B decays) 
and separately generated MC samples representing each individual search mode. The signal 
region for each search mode was then defined according to a cut on the output of the 
corresponding discriminator function. For each search mode, the value of this cut was 
selected in an unbiased way by minimizing the average expected MC Poisson upper limit V 
according to 

oo 

V = J2Piu,t)Bn{e) 

i=0 

where P{u, i) is the Poisson probability for finding i background events given an expectation 
of M, and Bri{e) is the 90% CL upper limit for the branching ratio if i events are found. The 
expected signal efficiency from the MC simulation at these optimal points ranged between 
24.8% and 37.9% for the various search modes, with expected backgrounds of between 0.0 
and 0.48 events. The efficiency is that for all B meson signal events, regardless of whether 
the decay occurred in the fiducial region of the detector, but does not take into account 
the branching ratios into fully-charged two-body final states for the vector mesons in the 
B ^ PV andB ^ VV modes. 



The signal efficiencies were determined from the SLD MC simulation, and thus are sub- 
ject to modeling uncertainties. The efficiency of the SLD tracking system was constrained 
by studying the track multiplicity distributions of inclusively tagged Z° — > bb events, which 
are identified with approximately 98% purity by the SLD [^. The kinematic distributions 
of tracks from such events are well constrained by measurements of B meson decay at 
the T(45') [[l^], as well as measurements of heavy-quark associated multiplicity at the Z^ 
pole |]TT|. The resulting comparison of the momentum dependence of the multiplicity be- 
tween inclusively tagged MC and data events indicated a deficit of ~ 5% in the tracking 
efficiency below 1.5 (0.8) GeV/c for the VXD2 (VXDS) data sample, leading to a reduction 
in the estimated signal mode efficiency of 6e/e ~ 1 — 2%. 

The possibility of longitudinal polarization of the vector mesons in the VV decay modes 
has been considered. A longitudinally polarized vector meson will decay with a cos^ 6h dis- 
tribution, with tracks from vector mesons decaying with small 6 tending to be reconstructed 
less efficiently, due to the relatively low momentum of the backward-going track, as well as 
the resulting angular proximity of the two decay tracks. The signal MC assumes 50% longi- 
tudinal polarization for the vector mesons from VV decays. Assuming a uniform probability 
distribution between 0% and 100% polarization, the resulting relative systematic error in 
the efficiency of the VV decay mode reconstruction is Ae/e = 1.5%. 

The momentum resolution at high momentum was studied by comparing the width of the 
reconstructed mass peak between data and MC for a sample of exclusively reconstructed 
D^ -^ K^-K^-K^ decays. To account for the somewhat larger width observed in data, 
the MC momentum distribution was smeared according to 1/p_l — > 1/p_l + Gaussian, for 
a Gaussian width of 0.002 (0.001) (GeV/c)-i for the VXD2 (VXDS) data sample. The 
resulting change in the MC mass width, for example for the B^ -^ p^n^ search mode, is 
from 146 to 184 MeV/c^ for the VXD2 data sample. The resulting efficiency loss varied 
between Se/e c^ 2 — 5%. Smearing of the radial and longitudinal track origin parameters, 
constrained by comparisons of r — and r — z impact parameter distributions between 
data and MC, yielded an additional efficiency loss of 6e/e ~ 2 — 4%. As a cross check, 
after the inclusion of the above corrections in the MC efficiency calculation, the number of 
reconstructed D^ -^ K^-k^-k^ decays is within S% of the MC expectation, well within the 
experimental uncertainty on the D^ —>■ K^ti^ti^ branching fraction and the D^ production 
rate. 

Application of the various search mode selection algorithms to the full 199S-8 SLD data 
sample yielded a total of four distinct candidate events {Ei - E^) which populated the signal 
regions of six separate search modes. The events observed (background expected) in each of 
these modes were as follows: event Ei for B^ —>■ p^p^ (O.Sl); events Ei, E2 for B^ ^ K p^ 
(0.49); events Ei, E2, E^ for 5° -^ K*°K*'^ {0.27); event E4 for fi° -^ (fK*\o.U); event E^ 
for B^ — > K* p°(0.S4); and events Ei, E3 for 5° -^ 'K* K*°(0.17). For the remaining search 
modes, no events were seen. 

Thirteen of fifteen MC events which passed the full selection criteria for at least one 
of the search modes had at least one identified track coming from a B meson decay, with 
B -^ D-K accounting for approximately one half of these. In four of the thirteen cases, 
reconstructed rest mass missing due to undetected charged or neutral particles was supplied 
by random fragmentation tracks. 

Figure 1 shows the expected signal and background distributions as a function of the 



discriminator output from the MC simulation of a typical mode - the B^ -^ p^p^ mode 
for the VXD-2 running period - assuming a branching fraction Bt{B^ -^ P^P^) = 10"*^. 
Figure 2 shows the relative rate between data and MC of the inclusion of background as 
the discriminator cut is relaxed, for the same sample. At a branching ratio of 10~^, a clear 
signal is expected, while backgrounds seem to be well modeled. 

The mode for which the observed signal was least likely to be accounted for by a statistical 
fluctuation in the expected background was B^ ^- K K*^. The Poisson likelihood of an 
expected background of 0.27 events fluctuating to three or more events is 0.27%, but depends 
strongly on the value of the expected background. A study for this mode similar to that of 
Fig. 2 yielded an additional 11 events, compared to a MC expectation of an additional 3 
events. Thus, for this mode there is reason to believe that the background is underestimated, 
and so, as for other modes, only an upper limit will be quoted. 

The branching ratio upper limits L are related to the statistical upper limits a on the 
number of observed events according to a = S ■ L; S = Nb ■ e, with A''^ the estimated 
number of applicable B meson decays, and e the estimated efficiency for reconstructing 
the given signal mode. The number of B^{B^) and B^ meson decays in the full SLD 
data sample is estimated from the measured SLD sample luminosity and known B meson 
production rates to be (1.02 ± 0.05) x 10^ and (0.27 ± 0.05) x 10^, respectively. It has been 
assumed that (21.7 ± 0.1)% of hadronic Z° decays involve primary b quarks, and of these, 
(39.7+^:^)% are 5° or B+ decays, and (10.5+^;^)% are 5° decays [||. The uncertainty in 
the signal mode efficiencies was conservatively estimated to be the total difference in the 
MC efficiency estimate with and without the extra momentum, tracking efficiency, and track 
origin parameter smearing. Including the additional VV mode polarization systematic error, 
as well as the MC statistical error of Ae/e ~ 2 — 5%, due to the limited size of the generated 
signal mode samples, the total modeling error was between 6 - 10% for all modes. 

Table 1 exhibits the number of candidate events, expected background, efficiency, sensi- 



tivity (S), and resulting 90% CL upper limits for both the Bayesian |T^ and Classical |T^ 
approaches for the 24 search modes. Each four-prong mode limit presented here, with the 
exception of B^ -^ (pK and B^ -^ 00, either improves upon the existing limit [^], or is the 
first available limit for the given mode. Two of the two-prong B^ modes (tt+tt^ and K~7t'^) 
are competitive with existing limits 0. Furthermore, a comparison of the probabilty distri- 
bution for the combination of the p^Tr"*" and K*^7[~^ modes with that implied by the DELPHI 
result BR(5+ -^ p^tt+, K*°tt+) = (1.7tJ;^ ± 0.2) x 10"^ yields only a 10% probability that 
the two measurements are consistent with the same central value. 

In conclusion, the excellent tracking capabilities of the SLD detector have enabled the 
SLD to establish a number of unique or competitive limits on the decay of B mesons to 
exclusive charmless final states. In particular, most of the four-prong quasi two-body lim- 
its presented here are the most stringent limits available. In addition, the SLD limits of 
BR(S+ -^ p^TT+) < 0.83 X 10"^ and BR{B+ -^ K*%+) < 1.19 x 10"^ (90% CL) rule out 
a DELPHI observation of the sum of these two modes [^ in favor of more stringent limits 
from CLEO |]. 
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TABLES 



Mode 


i 


a 


S (xlO-^) 


Bckd 


Data 


UL (B; xlO^) 


UL (C ; xlO^) 


B^- 


> vr+vr^ 


.338 


3.46 ± .31 


0.03 





0.69 


0.67 


B^- 


> K"7r+ 


.345 


3.53 ± .32 


0.14 





0.67 


0.66 


B^- 


yK+K- 


.341 


3.49 ± .28 


0.14 





0.67 


0.66 


B'.- 


> iT~^7r~ 


.379 


1.02 ±.16 


0.03 





2.35 


2.32 


B'.- 


>K-TT+ 


.335 


0.91 ±.15 


0.10 





2.62 


2.61 


B's- 


>K+K- 


.311 


0.84 ±.14 


0.20 





2.77 


2.83 


B+ - 


^p^7r+ 


.272 


2.78 ± .26 


0.34 





0.81 


0.83 


B+ - 


^p^K+ 


.264 


2.70 ± .24 


0.41 





0.83 


0.86 


B+ - 


^ K*\+ 


.285 


1.94 ±.17 


0.17 





1.21 


1.19 


B+ - 


^ K*^K+ 


.248 


1.69 ±.18 


0.17 





1.39 


1.38 


B+ - 


-> (f)1T'^ 


.301 


1.51 ±.12 


0.07 





1.59 


1.53 


B+- 


^(j)K+ 


.321 


1.61 ±.12 


0.14 





1.47 


1.44 


50- 


>p^p^ 


.263 


2.76 ± .24 


0.31 


1 


1.57 


1.36 


^0- 


> K*OpO 


.253 


1.76 ±.15 


0.49 


2 


3.30 


2.86 


B^- 


> K*^K*^ 


.304 


1.42 ±.12 


0.27 


3 


5.27 


4.69 


B^- 


>4>p^ 


.298 


1.53 ±.13 


0.14 





1.58 


1.56 


B^- 


>(j)K*^ 


.295 


1.01 ±.08 


0.14 


1 


4.34 


3.84 


B^- 


■^cfxt) 


.393 


0.74 ±.05 


0.00 





3.37 


3.21 


B'.- 


>p'p' 


.277 


0.77 ±.13 


0.27 





3.06 


3.20 


B^- 


> K*''p'^ 


.272 


0.50 ± .09 


0.34 


1 


8.52 


7.67 


B^- 


> K*^K*° 


.265 


0.33 ± .05 


0.17 


2 


18.21 


16.81 


B'.- 


>ct>p' 


.290 


0.39 ± .06 


0.07 





6.24 


6.17 


B'.- 


>(t)K*^ 


.265 


0.24 ± .04 


0.14 





10.02 


10.13 


B's- 


^4>4> 


.308 


0.21 ±.03 


0.00 





12.11 


11.83 



TABLE L Summary of efficiency (e), sensitivity (S), expected background, number of events 
in the signal region, Classical (C) and Bayesian (B) 90% CL Upper Limit for the 24 search modes. 
Note that the sensitivities (but not the efficiencies) take account of the branching fraction for (j) or 
K*^ into a fully-charged two body final state, where applicable. 
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FIGURES 
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Discriminator Output 

FIG. 1. Expected signal (shaded) and background for the B^ -^ p^p^ channel, as a function 
of the discriminator output. The plot shown is from the VXD-3 period Monte Carlo, assuming a 
branching fraction Br(B^ -^ P^P^) = 10~^. 
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Monte Carlo Events 

FIG. 2. Comparison between VXD-3 data and Monte Carlo of the rate of introduction of 
background into the B^ -^ p^p^ sample as the discriminator function cut is relaxed. The lines 
represent the upper and lower 90% CL limits under the assumption that the Monte Carlo accurately 
models the background. 
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